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Abstract 

A new coordinate system is defined for the Four-Body dynami- 
cal problem with general masses, having as its origin of coordinates 
the center of mass. The transformation from the inertial coordinate 
system involves a combination of a rotation to the principal axis of 
inertia, followed by three changes of scale leading the principal mo- 
ments of inertia to yield a body with three equal moments of inertia, 
and finally a second rotation that leaves unaltered the equal moments 
of inertia. These three transformations yield a mass-dependent rigid 
orthogonal tetrahedron of constant volume in the inertial coordinates. 
Each of those three linear transformations is a function of three co- 
ordinates that produce the nine degrees of freedom of the Four-Body 
problem, in a coordinate system with the center of mass as origin. 

The relation between the well known equilateral tetrahedron solu- 
tion of the gravitational Four-Body problem and the new coordinates 
is exhibited, and the plane case of central configurations with four 
different masses is computed numerically in these coordinates. 

Keywords: Four-Body Problem. New coordinates. 

PACS 45.50.Pk Celestial mechanics 95.10.Ce Celestial mechan- 
ics(including n-body problems) 
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1 Introduction 



The coordinate system introduced in this paper is a generahzation 
of the symmetric coordinate system of Piiia and Jimenez [1], [2], [3], 
that was defined for the Three-Body problem. Relative symmetric 
coordinates in the Three-Body problem were defined by Murnaghan 
[1] and Lemaitre [5] . More recently, Hsiang and coworkers have studied 
at least since 1995 the triangle geometry of this problem [6], with an 
important impact in the modern Three-Body problem reviewed by 
Chenciner [7J who posted in the web an important panorama on the 
subject, including the geometry of the so called shape sphere, that 
is almost the same coordinate system than ours for the case of three 
particles. Important contributions have been also made by Littlejohn 
and Reinsch [8] for the analysis of coordinate systems of three and 
four particles. 

The proposal of new coordinates, presented in this paper, have 
important points of contact with those works, although it sets itself 
apart from them, and simplifies their ideas in the case of four particles. 



The masses of the four bodies mi, m2, m-^ and m4 are generally dif- 
ferent, and we consider them ordered by the inequalities mi > m2 > 
1713 > rrii. 

We transform from the inertial referential, to the frame of principal 
axes of inertia by means of a three dimensional rotation G that is 
parameterized by three coordinates, such as the Euler angles. 

In addition to this rotation three more coordinates are introduced, 
as scale factors R2, R3, where Ri, R2 and R3 are three distances 
closely related to the three principal inertia moments through 



h=fi{Rl + Rl), l2 = fi{Rl + Rl), and h = fi{Rl + R^) , (1) 



The size of the scale factors is given in terms of the principal mo- 
ments of inertia by the equations 



2 The New Coordinates 



where is the mass 
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With the first rotation and the change of scale the resulting four 
body configuration has a moment of inertia tensor with the three 
principal moments of inertia equal. A second rotation G' does not 
change this property. 

The cartesian coordinates of the four particles, with the center of 
gravity at the origin, written in terms of the new coordinates are 

^ Xl 

Vi 



where G and G' arc two rotation matrices, each a function of three 
independent coordinates such as the Euler angles, and where the Oj, 
hj and Cj are twelve constants forming three linearly independent 4- 
vectors a, b and c, in the mass space, orthogonal to the mass 4-vector 
m = (mi, m2, ma, m4): 





aimi + a2m2 + a^m^ + a4m4 = , 
6imi + h2m2 + ^3^3 + 64m4 = , 
cimi + C2m2 + Cams + C4m4 = . 



(5) 



We introduce the following notation for the matrix 

/ mi \ 






M 



(6) 





m3 
m4 j 

In order to complete the definition of vectors a, b and c we assume 
aMb'^ = 0, bMc'^ = 0, cMa^ = 0, (7) 



V 





m2 





aM2b'^ = 0, hM^c^ = {), cM^a'^ = (8) 

which determine the directions of a, b and c in the 3-plane orthogonal 
to m, and we assume the normalizations 

aMa'^ = bMb^ = cMc^ = (9) 

that make vectors a, b and c adimensional. 

These three vectors are easily computed from the previous orthog- 
onality conditions. One has 

/I 1 1 1 \ 

a = , , , , (10) 

ymi-Xa m2-Xa ms-Xa m^-XaJ 
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/ 1 1 1 1 \ 

D = m , , , , (11) 

\mi — Xb 1712 — Xb — Xb m^ — Xb/ 

and 

/I 1 1 1 \ 
c = , , , , (12) 

\rai—Xc m2 — Xc m^ — Xc m^ — XcJ 

where i/a, Ub and t/c are normalization factors, and x^, Xb and x^ are 
the roots of the cubic equation 

— x^(mi+m2+m3+m4)+2x^(mirn2+rn2m3+?iT.3mi+m4mi+m4m2+m4m3) 

— 2>x{m2mzmi + m^m/^mi + mAmim2 + mim2m^) + ^mim2m2,rn4, = 0. 

(13) 

The symmetric nature of this equation is the consequence that this 
cubic polynomial is related to the derivative of the polynomial {y — 
■:^){v — :^)(v — ■:^){v — -i^)- The roots of this derivative are: I/xq, 
l/xb, l/a^c, and are located between the inverses of the masses. 

These quantities are defined in this form only for different masses. 
In that case we have the inequalities 

mi > Xa > 7712 > > ?7l3 > Xc > m4 , (14) 

that imply 

ai > 0, 02 < 0, 03 < 0, a4 < ; 

hi > 0, 62 > 0, 63 < 0, 64 < ; (15) 

ci > 0, C2 > 0, C3 > 0, C4 < . 

The column elements of the constant matrix 

ai 02 03 04 
E = I 61 62 63 64 I , (16) 

Ci C2 C3 C4 

are the coordinates of the four vertices of a rigid orthogonal tetrahe- 
dron. 

An orthogonal tetrahedron has the property that the perpendicular 
lines to the faces trough the four vertices intersect at the same point. 
Orthogonal tetrahedra were considered by Lagrange in 1773 [9j. Other 
old references on orthogonal tetrahedra are found in a paper by Court 
|10j . where he calls it orthocentric. Placing the four masses at the 
corresponding vertices, that intersection point is actually the center 
of mass of the four masses, and the moment of inertia tensor of the 
four particles has the same principal value in any direction. Equations 
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((Tj) and ([9]) imply that the inertia tensor of the rigid tetrahedron is 
proportional by a factor 2// to the unit matrix. 

To show these properties we consider a four vector linearly inde- 
pendent to the three 4-vectors a, b and c 



d=./^(l, 1,1,1) 
V m 



(17) 



where we use the notation m = mi + m2 + ms + rui for the total mass 
of the system. Then using definitions ([5]), d?]), ([9|), and (fT7j) we write 
them in terms of r = i/™ in the form 



1 



/ 



ai a2 as 
hi 62 ^3 

Cl C2 C3 

tj^ 



0-4: \ 

r J 



M 



( Ol 
V 04 



61 
64 



Cl r \ 

C2 r 

C3 r 

C4 r / 



/ 1 \ 

10 

10 

V 1 y 



(18) 

Since the inverse matrix from the left is equal to the inverse from the 
right, this equation transforms into 



/ ai 02 03 a4 \ I ^ 



( ai bi Cl r \ 

a2 62 C2 r hi 62 &3 ^4 

03 ^3 C3 r Cl C2 C3 C4 

\^ 04 64 C4 r y y r r r r 






V 





JL. 

1712 









JL. 

"13 





\ 





— y 



(19) 

Because this matrix equation is equal to its transposed; it just has ten 
independent equations. Four of them are 




The other six are 



aiGj + bihj + CiCj 



m 



(i = 1,2,3,4) 



(20) 



(21) 



From these basic equations it is easy to show that the position vector 
of one vertex is orthogonal to the three vectors between two vertices 
of the corresponding face (to the first vertex.) 

ai{aj - ttk) + hi{hj - bk) + Ci{cj - Cfc) = (i, j, k different) . (22) 

In addition, the distance between two vertices is given by 



1 



1 



+ {hi - hjY + (ci -CjY = ^^\— + — 



mi 



(23) 
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This is the condition to have a moment of inertia tensor with the same 
three principal moments of inertia. The six edges of the tetrahedron 
should be equal (proportional) to the square root of the right hand 
side of this equation. The volume of this tetrahedron is equal to 1/6. 

There are other remarkable geometrical properties of an orthogonal 
tetrahedron. The center of mass of each face is at the orthocenter 
where the three altitudes of the face intersect. This point is on the 
same straight line between the opposite vertex and the center of mass. 
In addition to the orthogonality of the three sets of two opposite edges 
of the tetrahedron, the two orthogonal edges are also orthogonal to 
the line joining the center of mass of the two edges. 

Let me do in this paragraph a technical digression that is specially 
relevant for engineering and physical minds. In formulating the ex- 
plicit expressions of the coordinates of the constant rigid tetrahedra E 
the origin for computing the G' rotation was arbitrarily chosen to be 
the one associated with the equilateral tetrahedra with four different 
masses, which implies a constant G' that was here selected equal to the 
unit matrix. This convention is introduced through equations ([8]) that 
are actually not necessary for the rest of the statements and proofs in 
this paper. Although there are of course other important coordinate 
systems to fix the origin for measuring the G' rotation; from these I 
prefer to choose one particle along one coordinate axis, the other three 
in a parallel plane to the parallel coordinate plane which does not in- 
clude the first particle; a second particle on an orthogonal coordinate 
plane that includes the first particle, and the other two particles on 
a line that is parallel to a coordinate axis and perpendicular to the 
coordinate plane of the first two particles. Another equally important 
referential for the origin of the rigid tetrahedron is associated with the 
grouping of the four particles in two sets of two particles. The center 
of mass of the two pairs, and the center of mass for the whole system 
are on a coordinate axis, and each of the two selected pairs of particles 
are placed on a line parallel to a coordinate axis. 

The previous definitions do not work in the important cases when 
two or more masses have exactly the same value. In those cases the 
tetrahedron is identified more easily from condition (j23p in terms of 
the masses. The selection of the origin for measure the rotation G' 
is now forced by the symmetry of the tetrahedron introduced by the 
mass equality. 

This rigid tetrahedron is the generalization of the rigid triangle of 
the Three-Body problem with the center of mass at the orthocenter 
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discussed previously in [TT]. 

I assume for simplicity that the potential energy is given by the 
Newton potential (the gravitational constant is equal to 1) 



V 



E 

Kj 



minij 



(24) 



although our results may be generalized for any potential which is 
a given power law of the relative distances between particles rij. It 
follows the relation between the interparticle distance and the new 
coordinates. The relative position between particles i and j is 



(25) 



The square of this vector is not a function of the first rotation G, but 
just of the scale matrix and the second rotation matrix 
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Xj 


— Xi \ 




( Rl 










Vj 


- Vi 


1- 
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V 




- Zi ) 




I 





Rs 





( flj 


- flj \ 






-k 






- Ci J 



( an 



^ij = ("-j ~ "-i bj - bi Cj - 
where A is the symmetric matrix 

= G' 



(26) 



/ All Ai2 ^13 
Ai2 A22 ^23 
\ ^13 ^23 ^33 



/ Rl 


V 



R2 









G 



/T 



(27) 



The six distances are thus functions of six components of matrix A or 
equivalently, are functions of the six independent coordinates in the 
scales Ri, and the rotation G'. 

We also compute the kinetic energy as a function of the new coor- 
dinates, which is given by 



2 



Ri^ - 4:{R2RiUJiVLi + R^Riuj2^2 + RiR2^M + 

( i?2 + -^3 



U=l 



\ 










Rl + R\ 
Rj + Rj 
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V 










, + Rt 







Rl + Rl 



(28) 



where uj = {uji,uj2-,i0'i) is the angular velocity vector of the first ro- 
tation G, and VL = (Oi,r22,^23) is the corresponding angular velocity 
vector of the second rotation G'. 



3 Equations of Motion 



The equations of motion follow from the Lagrange equations derived 
from the Lagrangian K — V as presented in any standard text on 
Mechanics [T2]. p^. 

However, the three coordinates related to the first rotation pro- 
duces Lagrange equations that imply, when the potential energy is a 
function only of the distances, conservation of the angular momentum 
vector in the inertial system 



G L 



(29) 



where L is the angular momentum in the principal moments of inertia 
frame 

dK 



duj 



{Rl + RlPi 
(Rl + Rl)u;2 
[Rl + Rl)uj^ 



2^1 



i?2^3^1 
R^RlQ2 
RlR2^i 



(30) 



This conservation leads to three first order equations forming, for 
this four body problem a generalization of the Euler equations valid 
for the rotation of a rigid body, namely 



/ /x(^i+i?i)a;i 
/x(i?i + RI)lo2 
\ fi{Rl + RDlos 



dt 



2^/22-^3^^! 

2fiR^RiQ,2 
2fiRiR2n3 



f fl{Rl - Rl)iJ2UJ3 + 2flRi{R2UJ2^3 - i?3C^3f^2) 

fi{Rl - Rl)u;3LJi + 2fiR2iR3i03ni - Riuin^) \ ■ (31) 

y H{Rl - Rl)uJlUJ2 + 2flR3{RiLOin2 - R20J2^l] 

The so called elimination of the nodes in the Three-Body problem 
, has a similar representation in this coordinates for the Four-Body 
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problem by means of the equation that equals the angular momentum 
vector in the principal moments of inertia frame to the rotation of a 
constant vector, which may be written in terms of two Euler angles 

/ + Rl)iOi \ ( R2R^^i \ ( 0\ 

III {Rl + Rl)uj2 - 2/^ R3Ri^2 = , (32) 

\ {Rj + RDcjs J \ RiR2^3 ) \l J 

where I is the magnitude of the conserved angular momentum. 

The Lagrangian equations of motion for the three scale coordinates 

are 

„ _ „ „ „ , „ . o o „o „o. dV 



f^-^Ri + 2n[R2L03n3 + R3UJ2n2] + nRi{u^ + a;| + + O^) - , 

(33) 

(f dV 

n-^R2 + 2n[R3UJini + Riuj3n3] + nR2iujl+(^j + nl + nl) = 

(34) 

and 

At^^3 + 2l^[RiL02^2 + R20Jl^l] + IJ.R^i'^l + + ^l) = ■ 

at 0R3 

(35) 

The three equations of motion for the three coordinates associated 
with the second rotation G' are written as an Euler equation similar 
to the one found for the first rotation, although the internal angular 
momentum is not conserved because of the presence of an internal 
torque 



di 



j ii{Rl + i?|)Oi - 21XR2R3UJ1 
fiiRl + Rl)n2 - 2iiR^RiU2 
\ n{Rl + Rl)n^ - 2fiRiR2UJ3 




( n{Rl - i2f)f^2J^3 - 2nRi{R2U2^3 - R3i^3^2) \ 

fi{Rl - Rl)n3ni - 2iiR2{R3UJ3^i - RiuJiQ. 



3 J 



(36) 

\ IJi{Rl - Rl)VLiVL2 - 2nR3{RiLUin2 - R2^2^i) j 

where Ki,K2-, K3 are the components of the internal torque K which 
is expressed in terms of the derivatives of the potential energy with 
respect to the three independent coordinates qj in the rotation G' 
and the three vectors Cj that appear in the expression of the angular 
velocity in terms of the same coordinates 

3 

= ^c,g, (37) 
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where the vectors Cj are generally functions of the coordinates qj. 
The internal torque is determined by the equations 



K c. 



dV 



(38) 



There is one more constant of motion, namely the total energy 

3 



E = V+K = 



J2R^ - 4(i?2i?3'^lf^l + R?.RlOJ2^2 + RlR2^Z^z) + 



i=l 



Rl + Rl 

Rl + Rl 

Rl + Rl ) 

Rl + Rl 

Rl + Rl 

Rl + Rl 



iO + 



(39) 



4 The plane problem 



The case with the four particles in a constant plane is an important 
and old subject [15j. Our coordinates are now adapted to that case. 
The third component of the cartesian coordinates of the four particles 
are zero. The modification of our coordinates (4) for this case is given 
by two changes: the first rotation by just one angle in the plane of 
motion; and the scale associated with the third coordinate is zero, 
namely 

X2 X3 X4 



/ 



Xl 

2/1 




2/2 





2/3 




2/4 






^ ai 
hi 

V ^1 



This equation simplifies to 



a2 
h2 

C2 



03 
63 
C3 



a4 
64 
C4 



(40) 



Xl X2 X3 X4 

2/1 2/2 2/3 2/4 
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/ cosV' -simp \ [ Ri \ ,t ( 1 
sinV' cos^ i?, ' 

V / \ / \ Ci C2 C3 C4 / 

(41) 

in terms of six degrees of freedom. 

We need three independent coordinates (for example three Eu- 
ler angles) in G' for the two independent vectors in four dimensions 
expressed in the base of the three constant vectors a, b, and c, or- 
thogonal to the mass vector. 

Let me do in this paragraph a technical digression that is specially 
interesting for engineering or physical minds. We must formulate in a 
mathematical language the conditions for a plane solution. The most 
usual way to do this is to equal to zero the Cayley-Menger determi- 
nant which has entries equal to 1, 0, and the squares of the distances 
between particles. Although Dziobek [15] considered this approach of 
paramount importance, however he introduced equivalent conditions 
that have been promoted by many years by A. Albouy and cowork- 
ers (see [16] and references therein,) which consists in using the four 
directed areas of the triangles formed by the particles. 

The four (twice) directed areas are written in terms of the cartesian 
coordinates as 





1 


1 


1 
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1 


Si = 








, 5*2 — 


Xl 


X4 


x^ 




2/2 


2/3 


2/4 




2/1 


2/4 


2/3 
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1 


1 




1 


1 


1 


S3 = 


Xi 


X2 




, 5*4 = 


Xl 


X?, 


X2 




Vi 


2/2 


2/4 




2/1 


2/3 


2/2 



that are the four signed 3x3 minors formed from the matrix 

/ 1 1 1 1 \ 

Xl X2 X2, X4 (43) 

V 2/1 2/2 2/3 2/4 / 

Addition to the previous matrix of a row equal to any of its three 
rows produces a square matrix with determinant zero, that implies 
that the necessary and sufficient conditions to have a constant plane 
tetrahedron are 

E5, = 0, (44) 

i=l 
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and 



4 4 

^ SiXi = , ^ 



. (45) 
The two last equations summarized by the zero vector condition 

4 

Y.Siri = 0. (46) 

i=l 

An expression for the three directed areas in terms of the previous 
coordinates follows 



/ Si \ 

V ^4 J 



CME^G' 




(47) 




where C is a constant with units of area over mass. Substitution 
of equations (j4ip and (j47p in equations (j44p or (j46p one obtains an 
identity, independent of coordinates i?2, "0 c'-iid the rotation angle 
of G' around the unit vector 



(48) 



Given the four masses, the four directed areas of the four particles 
are functions of this unit vector direction only, up to a multiplicative 
constant C depending on the choice of physical units. These explicit 
expressions should make clear Albouy's |16j affine formulation of the 
plane condition. Another form of this constant plane condition is also 
published in reference [20j. 

In the plane case the angular momentum has a constant direction 
orthogonal to the plane and of magnitude 



dK 

dtp 



(49) 



The kinetic energy becomes 

2 



K 



2 



.i=l 



12 



( Rl 
ii? 








\ Rl + Rl 



(50) 



Substitution of polar coordinates for the Ri and R2 coordinates 

Ri = Rcose, R2 = Rsme (51) 

and writing the kinetic energy in terms of the angular momentum 
constant of motion instead of the ip velocity lead us to 



2 



if + R^ (e'^ + cos2(26l) + sin^ 6 + ^1 cos^ 9 
Energy conservation is thus expressed as 



+ 



p2 

2/ii?2 • 
(52) 



E 



R^ + R^ ((9^ + n| cos2(26l) + sin^ 6 + ^1 cos^ 6* 



p2 

(53) 



where 1^ represents the potential energy. 



5 Central configurations 

In this section we begin with the approach by Dziobek [15]. See refer- 
ence [T7j for a contemporary approach. The Four-Body central config- 
urations are determined as critical points of the potential energy with 
a fixed total inertia moment that in three dimensional space lead to 

mirrij/rfj = arriimj . (54) 

The left hand side of this equation is the derivative of the potential 
energy with respect to rfj . The right hand side is the derivative of the 
inertia moment with respect to rfj multiplied by an unknown constant 
a that includes the constant total mass, contained in the expression 
for the total inertia moment, and the gravity constant in the potential 
function. Of course, this equation simplifies to 

1/4 = a. (55) 

According to equation (j55p . the only Four-Body three dimensional 
central configuration results just if the six distances are the same, 
giving an equilateral tetrahedron. For an equilateral tetrahedron one 
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particular coordinate system is given placing its vertices on alternat- 
ing corners of a cube having the six faces normal to the coordinate 
axis. Then the center of mass is computed and the origin of coordi- 
nates translated to this position. After that the tensor of inertia is 
determined and the scale factors associated to the principal moments 
of inertia are in the ratios 

Rl/B^/Rl = xa/xi/xc . (56) 

This fact is deduced because the obey the same characteristic equa- 
tion (113p . Now we show that the equilateral tetrahedron has been 
selected as origin for measuring the G' rotation. A tetrahedron with 
positions 



\JXal\l 




I I 6i 62 &3 ^4 I (57) 

\ v^T^ 

is an equilateral tetrahedron (namely the G' rotation is the unit matrix 
for the equilateral case). The proof is obtained from this equation by 
direct computation of the length of the edges and use of equations 
(j20|) and (pT]) . The six edges of that tetrahedron are equal to 

The equilateral tetrahedron gives a well known solution in which 
the masses move on straight lines. This is a theorem by Laplace |18] 
with an old history. 

The non collinear plane central configurations are characterized 
in our coordinates by constant values of the G' matrix and of the 
coordinate Q associated to the constant value of the ratio -Ri/i?2- For 
these cases the angular velocity vector is the null vector, the angular 
velocity Q is also zero and the equations for conservation of moment 
and energy, (H9]l and (i53|) respectively, become 



= fiipR^ . (58) 



and 



These equations are identical to similar equations obtained for the 
Euler and Lagrange central configurations of the Three-Body problem 
[21j . They are formally the same as the equations for the conies in the 
Two-Body problem in terms of the radius R and the true anomaly "0- 
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The constant values of the G' matrix and angle refered to above 
are not arbitrary but they are determined by three independent quan- 
tities as discussed in the following. 

The plane solutions with zero volume but finite area are obtained 
taking in account that the variational equation (j54p is modified adding 
the restriction of plane motion. This condition is obtained by Dionzek 
|15j from the derivative of the Cayley-Menzer determinant with re- 



spect to rfj that he found to be proportional to the product of the 
directed areas SjS' 



3- 



It follows that the solution is given in terms of parameters A and 



a 



rjk = ^ + ^^i^k , (60) 

where Aj = Sj/rrij are weighted areas, quotient of the directed area 
divided by the corresponding mass. This equation was presented by 
Dziobek [15]. A proof was published by Moeckel [19], and using a 
different approach to the same problem, deduced by Albouy [T6]. A 
new proof of the equation was obtained in a different approach by 
Piiia and Lonngi |20j . 

It follows from (147p that in a plane solution the weighted directed 
areas are expressed as 



(A, \ 

A2 

Az 

V ^4 / 



/ 



\ 



(61) 



The weighted directed areas are up to a normalization factor equal to 
the third rotated coordinate of the rigid tetrahedra. 
The weighted directed areas obey the condition 



Y^Ajnij = 



(62) 



expressing the fact that the sum of the directed areas is zero, equation 



Since the lengths and masses are defined up to arbitrary units, I 
assume with no loss of generality that the parameter a equals unity. 



' jk 



1 + XAjAk 



(i / k). 



(63) 



This equation has been considered giving particular values of the four 
masses and computing the six distances by solving it under re- 
striction (j44|) . The weighted directed areas Aj are functions of the 
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2 

Figure 1: Two different sets of simultaneous positions of four particles with 
different masses following elliptic trajectories in a central convex plane con- 
figuration. The isolated point is at the center of mass at the common focus 
of the four ellipses. The eccentricity of the four ellipses is e = 0.72. 

distances and the masses only. Some examples of such approach are 
|16j . |22j . [23]. According to D. Saari [T7] the problem with this per- 
spective is difficult to manipulate, but we found it to be perfectly 
feasible as follows below. 

In the paper by Piha and Lonngi [20] a different point of view was 
adopted, namely that the directed weighted areas (that are defined 
with a simple functional dependence with respect to the masses,) are 
known as four given constants. The previous equation then gives the 
distances as functions of the unknown parameter A. Trough them, the 
areas of the four triangles become functions of A, that should obey 
the necessary restriction (I44p . to verify that one has a plane solution. 
This restriction allows in many cases to determine the value of A and 
hence the values of the six distances and the four masses. This is an 
implicit way to deduce planar central configurations with four masses. 

From the distances and masses one determines the positions of the 
four particles in the plane frame of principal moments of inertia and 
the principal moments of inertia are also computed. This allows to 
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know eight components of the rotated rigid tetrahedron EG' and the 
remaining coordinates are known from the four given weighted areas 
constants Aj according to equation ([6T]) 

Following this method we computed the necessary data to plot the 
trajectories of the four particles with different masses represented in 
Fig. 1. We started from the four constants 

Ai = 6, A2 = -15, A3 = 4, ^4 = -3 

The constant plane conditions give us the value 

A = -0.01268487093192263... 

from which the distances are 

r23 = 1.12863753386515... 
rai = 0.933745641175193... 
ri2 = 0.953868245971217... 
Ui = 1.32572435746881... 
r42 = 0.828080639336103... 
r43 = 0.775802698722361... 

and the corresponding masses 

mi = 0.428260218865972... 
ma = 0.355184464717379... 
m3 = 0.261905866491155... 
m4 = 0.113826160081235... 

This information is sufficient to compute an initial central config- 
uration and from it the four constant coordinates determining G' and 
6. The elliptic choreography (with no symmetry) is obtained simply 
by writing coordinate i? as a function of the real anomaly ip. The 
latus rectum, the eccentricity and the initial value of il^ may all be 
selected arbitrarily. 
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